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Abstract 

We microscopically evaluate the excitation spectrum of the ^He impurity 
in liquid ^He at T=0 and compare it with the experimental curve at the 
equilibrium density. The adopted correlated basis perturbative scheme in- 
cludes up to two independent phonon, intermediate correlated states and the 
correlation operator is built up with two- and three-body correlation func- 
tions. The experimental spectrum is well described by the theory along all 
the available momentum range. A marked deviation from the simple Landau- 
Pomeranchuck, quadratic behavior is found and the momentum dependent 
effective mass of the impurity increases of ~ 50 % at g ~ 1.7 A^^ respect to 

its g = value. No signature of roton-like structures is found. 
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Typeset using REVT^ 



Both experimentahsts and theoreticians have devoted a great deal of efforts to measure 
and explain the characteristics of one ^He impurity in atomic liquid "^He. From the exper- 
imental point of view it is well known how the impurity chemical potential, yUs, behaves 
with the temperature and the pre J, it. eflective ma==, and qua^iparticle excita- 
tion spectrum, e(g)lii. In ref.i the authors found a sizeable deviation from the quadratic 
Landau-Pomeranchuk (LP) spectrumH, eLp{q) = h'^q^/2ml, in low concentrations ^He-^He 
mixtures. e{q) was parametrized in a modified LP (LPM) form as: 

h^q'^ 1 

eLPAiiq) = TT-T 2 ■ (1) 

The estimated values of the LPM parameters, at P = L6 bar and xs ~ 0.05 (a:3=^He 
concentration), are rrig ~ 2.3 and 7 ~ 0.13 A^. 

Microscopic calculations, done in the framework of the correlated basis function (CBF) 
perturbation theory!, have been able to give good estimates of /is and ml at the T=0 
^He equilibrium density, peq = 0.02185 A~^. Lately, a diffusion Monte Carlo approach 
has provided similar results!. There are also theoretical indications of a deviation of the 
spectrum from the LP formiH. The presence of a possible roton-like structure in e{q) near 
the crossing with the ^He phonon-roton spectrum was supposed in ref.lHl but not confirmed in 
refs.i'i. However, in the last work an excitation spectrum quite higher than the experimental 
one was found. 

Here we will employ the CBF machinery of ref.i (hereafter denoted as I) to compute, 
in a microscopic way, the whole impurity spectrum. The CBF basis used in I consisted in 
correlated n-phonon states 

^q;qi..q„ = P3(q - Ql - .. - qn)p4 (Qi) • •P4 (Qn) ^0 (2) 

where P4(k) = J2i=i,N4^''^'^' is the "^He density fluctuation operator and P3(k) = e*^""^ de- 
scribes the excitation of the impurity. The basis states were then properly normalized. 
\E'o = \E'o(3,A^4) is the ground state wave function of '^He atoms plus one ^He impurity 
in a volume Q, taken in the Q 00 limit, at constant ^He density, = N^/Q. 



A realistic choice for \l/o(3, is made by applying an extended Jastrow-Feenberg cor- 
relation operatoilll to the non interacting g.s. wave function 

*o(3, N^) = F2(3, N^) Fs{3, N^) $o(3, N^) . (3) 

F2^3 are N-body correlation operators including explicit two- and three-body dynamical cor- 
relations. F2 is written as a product of two body Jastrow, ^He-^He and ^He-^He correlation 
functions, 

^2(3,iV4)= n /^'''H^3.) n f'''\rim), (4) 

i=l,N4 m>l=l,N4 

and F3 is given by the correspondent product of triplet correlations, f^"''^'"'\va,rp,Vj). 

The correlation functions are variationally obtained by minimizing the g.s. energy of the 
system, Eo. The procedure is outlined in I, where a parametrized form for the triplet correla- 
tions was used and the Jastrow factors were obtained by the Euler equations, SEo/Sf^"'^^ = 0. 
The equations were solved within the HNC framework and the scaling approximation for 
the elementary diagramjHI. The Aziz interatomic potentiafl was used in the minimization 
process. 

The perturbative calculation of I included one and two independent phonon (OP and 
TIP, respectively) states and all the diagrams corresponding to successive rescatterings of the 
one phonon states (OPR). This contribution was obtained by solving a Dyson-like equation 
in the correlated basis. While the correlation factors are intended to care for the short range 
modifications of the ground state wave function due to the strongly repulsive interatomic 
potential, the basic physical effect induced by the perturbative corrections may be traced 
back to the backfiow around both the impurity and the ^He atoms. The CBF analysis 
provided fi^iCBF) = -2.62 K (vs. j^^iexpt) = -2.79 K) and ml{CBF) = 2.2 mg at 
equilibrium density. The chemical potential was obtained with the Lennard- Jones potential 
and some improvement may be expected by the Aziz interaction. 

In order to construct the CBF perturbative series, we write e(g) = eo(g) + Ae(g), with 
eo{q) = hq'^/2m3 and 
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Ae(g) ~ Aeop(g) + Aeripiq) + A(^oPR{q) ■ (5) 

The different terms in represent contributions from the corresponding intermediate 
states. The n-phonon states have been Schmidt-orthogonahzed to states with a lower num- 
ber of phonons. For instance, the actual OP state reads as: 



„. „ \ _ l^q;qi/ I ^q/\^ql ^q;qi/ ( p\ 



l^q;qi/ 1 


^q)(^ql 


^q;qi/ 


\^q;qi 


IvT/ \l/2 

1 ^q;qi/ 



The two-phonon state, \l'q;q^q2, has been orthogonalized in a similar way to \E'q, ^q;qi+q2 
and ^E'q;qi 2 • The orthogonalization is a necessary step in fastening the convergence of the 
series as the non orthogonalized states have large mutual overlaps 

The non diagonal matrix elements (MEs) of the hamiltonian, H, (we recall that we use 
the Aziz potential) are evaluated by assuming that the two- and three-body correlations 
are solutions of the corresponding Euler equations. This is not strictly true for the triplet 
correlations but the corrections are expected to be small. With this assumption, it is easily 
verified that 

{q\H\q- qi) = -[N,Siq^)]-'/'^q ■ q^S^iq,) , (7) 

2m3 

where S{qi) and 5*3(^1) are the ^He and impurity static structure functions. 

In general, MEs involving n — 1 phonon states, are expressed in terms of the n-body 
structure functions 



and 



.(n). „ ^ (^o|P4(qi)--P4(qn-l)P3(qn)|^o) 



with q„ = qi + .. + q„_i. 

The diagonal MEs have a particularly simple form: 



(q; qi-qn|q; qi..q„) = N2S{q,)..S{q^) , (10) 
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and 



(q; qi-qn|^^|q; qi--qn) = + eo(g) + ^F^qi) (n) 

i=l,n 

with = (\l/o|-f:r|\l/o)/(^o|^o) and wpiqi) = h'^qf / 217148 ^{qi) is the Feynman ^He excitation 
spectrumO. 

The OP and TIP perturbative diagrams contributing to Ae(g) are shown in Fig. (5) of 
I, where only their q = derivative was computed, since the paper was concerned with just 
the calculation of the effective mass at g = 0. Here we extend the formalism to finite q. We 
use Brillouin-Wigner perturbation theory, so the correction itself depends on e{q) and the 
series must be summed self-consistently. For instance, the OP contribution is given by 

\{ci\H-Eo-eiq)\ci;q,)\' 



Aeop(g) = 



qi 



e(g) -eo(lq-qil) - ^^(^1' 



The expressions of the other diagrams are quite lengthy and will not be reported here. 
However, some comments are in order. They involve the two- and three-body structure 
functions, i.e. the Fourier transforms of the two- and three-body distribution functions, 
and (^'^^^(ri, r2, rs). g^"^^ is a direct output of the HNC/Euler theory and, in pure 
"^He, results to be very close to its experimental measure. To evaluate g^^^ is more involved 
and usually one has to resort to some approximations. The mostly common used are the 
convolution approximation (CA) and the Kirkwood superposition approximation (KSA)0. 
The CA correctly accounts for the sequential relation between g^^^ and g^'^^ and factor- 
izes in momentum space, 5'^j^(qi, q2, qs) = S{qi)S{q2)S{q3); the SA factorizes in r-space, 
(7j^5^(ri, r2, ra) = g'''^\ri2)g^^\ri3)g'^^\r23), and describes adequately the short range region. 
The momentum space factorization property makes the use of the CA particularly suitable 
for our perturbative study. 

The sensitivity of the calculation to the approximation for g^^^ clearly shows up in the 
CBF-TIP evaluation of the ^He excitation spectrum, u;(g). Figure 1 compares the Feynman 
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spectrum and those obtained within the CA and KSA with the experimental data. The 
phonon hnear dispersion at low-g is well reproduced by both Up^q) and ujca{.(i)i whereas 
^ksa{<i) fails to give the correct behavior. As it is well known, the remaining part of the 
spectrum is severely overestimated by upiq)', both CA and KSA give a reasonable description 
of the maxon region but KSA is closer to the experiments at the roton minimum, just in 
view of its better description of the short-range regime. An overall good agreement with 
the experimental curve was obtained in ref.0 where backflow correlations were added to the 
CBF states. 

Figure 2 shows e(g) in CA and KSA, along with the data from ref.i. The curves do not 
include the OPR contribution. At this level, the effective masses are m^{CA) = 1.6 
and mT^{KSA) = 2.1 and, again, KSA is closer to the experimental spectrum at large 
momenta. The curve labelled as CAl is obtained in CA, but using the experimental ^He 
spectrum in the energy denominators. Diagram (5.e) of I, that gives the two-phonon cor- 
rection to uj{q), has not been included as its effect is mostly taken into account by the use 
of ujexptiq)- KSA and CAl are close at large g- values, pointing to a good description of the 
^He roton as a key ingredient for a correct approach to the large q sector. We will follow 
the CAl method for the remaining of the work. 

Figure 3 gives the CAl impurity spectrum and the experimental ^He and ^He curves. The 
OPR terms are included and the LP and LPM fits to e^xptiq) are shown. Since the branch 
of the dynamical response due to the excitations of the low concentration ^He component in 
the Helium mixtures overlaps the collective ^He excitation at g > 1.7 A~^illi, Cexptiq) is not 
known in that region. A roton-like behavior was supposed in ref.i. This structure was not 
confirmed by the variational MonteCarlo (VMC) calculation of ref.i, which employed shadow 
wave functions in conjuction with a Jastrow correlation factor of the McMillan type. The 
VMC data at equilibrium density are given in the figure: they overestimate the experiment 
and have an effective mass of Tri^iV MC) ~ 1.7 m^. 

The shadow wave function of ref.i takes into account backfiow effects. Actually, in several 
papers was pointed out that second order perturbative expansion with OP states introduces 
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backflow correlations into the wave functionB'ElO. We find m\{OP) = 1.8 ms, in good 
agreement with the VMC outcome. An analogous CBF treatment by Saarelalll gave similar 
results (mg ~ 1.9 777,3) and a spectrum close to the LP form. More complicated momentum 
dependent correlations are generated by TIP and OPR diagrams, playing a relevant role in 
the CBF approach and giving m*^{CBF) = 2.1 m^. 

The total CBF impurity spectrum is very close to eexpt{(l) up to its merging into the 
^He dispersion relation. For the 7 parameter in the LPM parametrization, the theory gives 
'-){CBF) ~ 0.19 A?. If the spectrum is parametrized in terms of a momentum dependent 
effective mass, e(g) = ti^q^ /2mj^{q)., then we find 7773(g = 1.7 A^^) = 3.2 7713, with an increase 
of ~ 50% respect to the q = value. 

Beyond g ~ 1.9 the energy denominators vanish for some momentum values and the 
series cannot be summed anymore. This is due to the fact that the impurity quasiparticle 
is no longer an excitation with a well defined energy, since it can decay into ^He excitations 
and acquire a finite lifetime, r. A finite r-value reflects a non zero imaginary part of the 
■^He complex optical potential (or the on-shell self-energ y), W{q) = $5S(g,e(g))0. Figure 3 
shows W{q) as computed with only OP intermediate states, 

Wop{q) = vrE l(ql^ - ^0 - e(g)|q; qi)r5(e(g) - eo(|q - qi|) - w{q,)) , (13) 
qi 

where the LPM impurity spectrum and the experimental ^He dispersion have been used 
(notice that W{q) is amplified by a factor 4 in the figure). The OP optical potential is close 
to the one found in ref.lii. A numerical extrapolation of the computed ecBpiq) into the 
roton region does not show any evidences of a ^He roton-like structure. 

In conclusion, we find that CBF perturbative theory is able to give a quantitative de- 
scription of the ^He impurity excitation spectrum in liquid ^He at equilibrium density. The 
intermediate correlated states must consider at least two independent phonon states and 
one phonon state rescattering is found to play a non marginal role at large momenta. It is 
plausible that in a richer basis, including for instance explicit backflow correlations, a lower 
order expansion might be sufficient. However, the more complicated structure of the matrix 
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elements could cause additional uncertainties in their evaluation, at least in the framework 
of the cluster expansion approach. The development of a MonteCarlo based alghoritm for 
the computation of non diagonal matrix elements would probably be the correct answer. 
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FIGURES 

FIG. 1. ^He excitation spectrum at equilibrium density. Stars are the experimental data. 
Energies in K and momenta in A~^. 

FIG. 2. ^He single particle energies in CA, KSA and CAl without phonon rescattering. Stars 
are the experimental data. Units as in fig.(l). 

FIG. 3. CBF/CAl (triangles), LP and LPM ^He single particle energies. Full circles are 
the VMC data. Stars and circles are the impurity and ^He experimental data, respectively. Black 
triangles are extrapolated CBF/CAl values (see text). Black diamonds give the impurity imaginary 
optical potential (in K). Units as in fig.(l). 
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